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Abstract. Orthorhombic SrIrOs subjected to strain show tunable transport properties. With 
underlying symmetry remaining invariant, these properties are associated with IrOe octahedral 
tilting. Adopting to first-principles methods, the effects of crystal field, spin-orbit coupling, 
and Coulomb correlations, on comparable interaction length scales, are discussed. While tilting 
induces a t 2 g — eg crystal-field splitting and band narrowing, spin-orbit coupling induces a partial 
splitting of the Jeff bands rendering SrIrOa a semi-metallic ground state. The SOC enhanced 
hybridization of Ir-O orbitals, serve as a explanation to why the critical Hubbard correlation 
strength increases with increasing SOC strength in SrIrOs to induce an insulating phase. 
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1. Introduction 


In the Ruddlesden-Popper series Sr„+ilr„03„+i, n = 1 (Sr^IrO^) fl |-3 and n = 2 (81311207) 
[M 3 are insulators, while n = 00, (SrIrOs) is semi-metallic All systems in this series, 

chemically associate Ir with a formal valence of -|- 4 , due to its octahedral coordination (IrOe) with 
the neighboring ligands. The insulating ground state of Sr2lr04 and Sr3lr207, with Ir in its bdf 
electronic configuration, is accounted by means of the spin-orbit coupling (SOC) driven Je// model. 
According to the model, the t2g bands in Sr2lr04 and Sr3lr207 are split into Je// = | and Je// 
= ^ multiplets The Jeff = I states accommodate four electrons and therefore, are 

completely filled. The remaining electron is occupied into the Jeff = 5 doublet state. However, by 
virtue of their narrow bandwidth Coulomb correlations set in, which splits the Je /f = \ doublets 
into two energetically distinct Hubbard bands, thereby, leading to an insulating ground state [ii- 


Interestingly, despite having a octahedral coordination with a formal Ir^’*' state, SrIrOs has 
a semi-metallic ground state. Fundamentally, this may be attributed to the underlying crystal 
structure, itself. Both Sr2lr04 and Sr3lr207 crystallizes in tetragonal structure in which the IrOe 
motifs are stacked along the crystallographic c—axis, separated by Sr ions. The IrOe octahedra 
are distorted and tilted about the c-axis by an angle 11 ° in Sr2lr04 and 12 ° in Sr3lr207[3 24 , 
Since there is no direct linking of the these motifs along the c—axis, both Sr2lr04 and Sr3lr207 
display a quasi-two dimensional structure, commonly referred as single layered and bi-layered 
iridates , respe ctively. On the other hand, Srlr03 crystallizes in both monoclinic and orthorhombic 
phases | 2 bl 428 j . with the former being the ambient phase. The orthorhombic structure is stabilized 
when synthesis conditions are adopted to high pressure and/or when grown in thin film form on 
appropriate substrates [HI. [2^ 34 1. In both monoclinic and orthorhombic phases, the IrOe motifs of 
Srlr03 are interlinked along all three crystallographic axes, thereby imparting a three dimensionality 
to its structure. Thus, unlike Sr2lr04 and Sr3lr207, there exists an additional Ir-O-Ir chemical 
bonding in SrIrOs. Due to the additional structural constraint of the ions imposed by the c—axial 
bonding in Srlr03, the IrOe octahedral tilt is enhanced to 16 ° [l^, partly accounting the steric 
displacement of the Sr ions. 

Octahedral distortions in perovskites are generally associated with size mismatch of the 
constituents and the nature of chemical bonding. Such distortions are manifested as tilting and/or 
elongation of the octahedra. Elongation usually results due to the dominant breathing modes which 
originate due to charge disappropriation in the system, as like in BaBi03 where Bi exists in both -|-3 
and - 1-5 valence states. However, multiple valence states associated with Ir in Srlr03 seems less likely, 
and thus octahedral tilting modes appear more dominant. Therefore, the energetics associated with 
the structure of Srlr03 can be associated with the tilting of the IrOe octahedra. Identifying the 
IrOe octahedral tilt angle ((/) as a physical parameter, we study the electronic structure properties 
of SrIrOs using first-principles density functional theory based methods. We find that octahedral 
tilting not only induce a t2g - Cg crystal field splitting but also leads to significant band narrowing. 
Further, one finds that inclusion of SOC in the crystal Hamiltonian effectively splits the low energy 
antibonding states into Jeff states which are strongly hybridized. The calculations predict a semi- 
metallic ground state, which are consistent with the experimental observations. Also, following 
the empirical rigid band model of the alloy theory, it is anticipated that the transport properties 
of Srlr03 would sensitively depend on the concentration of charge carriers, as well. Furthermore, 
adopting to Hubbard Hamiltonian, our calculations also show that Coulomb correlations in SrIrOs 
have marginal effects on its low temperature electronic structure properties. 
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2. Computational details 

Optimization of the structural parameters and the electronic structure properties of SrIrOa are 
carried out using the full potential linearized augmented plane-wave (FP-LAPW) method as 
implemented in the Wien2k code [s^. The LAPW sphere radii for Sr, Ir and O were chosen 
as 2.26, 2.10 and 1.70 a.u., respectively. The lattice parameters were adopted to the experimental 
values, with a= 5.591 A, b =5.561 A and c = 7.882 A [^. The internal ionic coordinates of the 
Sr and O ions, which are not fixed by the symmetry, are relaxed using the force minimization 
technique. Well-converged basis sets were ensured, by choosing the Wien2k parameters RKmax= 

8 ; Gmax = 24 and Imax = 10. Additional local orbitals were also used to account for the semi-core 
Ir 5p states. The exchange correlation potential to the crystal Hamiltonian was considered in the 
generalized gradient approximation (GGA) as prescribed by Perdew, Burke and Ernzerhof [^ . 
The Brillouin zone fc—mesh sampling for total energy convergence was modeled with 384 fc-points 
in its irreducible part. The core states were treated relativistically, while SOG was included for the 
valence states through the second variational step (3^ . 

3. Results and discussions 
3.1. Structure 
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Bond distance (A) Bond angles (°) 


Ir-O(l) Ir-0(2) 
Expt 2.0215 2.0281 
Theory 2.0251 2.0303 


Ir-0(1)-Ir Ir-0(2)-Ir 

154.205 152.491 

153.346 152.399 


Table 1. The GGA optimized structural parameters of SrlrOa - the ionic coordinates, selected 
neighboring bond distances and bond angles, compared with the experiments l28ll . 

To estimate the internal coordinates of Sr and O ions, force optimization was performed with 
its tolerance set to 1 mRy/a.u. The resulting structural parameters are found in good agreement 
with the experimental data [2^. The theoretical and experimental values are compared in Table [T] 

The total energy minimization were also carried out by varying (j) to ascertain whether SOC has 
any prominent role in determining the structural properties of SrIrOs. The values of </> corresponding 
to these non-equilibrium structures, were estimated by using the relation cosOi = and 

Here, 9i and O 2 are the Ir-0(1)-Ir and Ir-0(2)-Ir bond angles with 0(1) and 
0(2) representing the out-of-plane and in-plane O ions of the orthorhombic unit-cell. The average 
of (pi and p 2 (= </>) was used as the parameter to minimize the total energy. The variation in 
the total energy of SrIrOa with respect to tilt angle, E{(p), is shown in FiglTJ Compared to the 
experimental value of 15.8°, we find that the scalar relativistic GGA calculations overestimate the 
equilibrium (p as 16.5°, while the GGA+SOC calculations underestimate it to 15.2°. 

The origin of octahedral tilting in perovskites is mainly attributed to steric and/or electronic 
interactions. In the steric models, tilting is facilitated due to the misfit size of the cation. Following 
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Figure 1. Color online ; The relative change in the total energy, with respect to the equilibrium, 
as a function of 0 in SrIrOs, calculated in the scalar relativistic GGA scheme ( blue line) and 
with SOG (red line) included in the Hamiltonian. 


the GGA and GGA+SOC total energy calculations, the difference in the estimated equilibrium value 
of (j) differs more than 1°. Also, from FigH] it follows that the stiffness of the tilting, a quantity 
which is proportional to the second derivative of E((j)) curve at equilibrium, increases substantially 
with SOG. The estimated increase with respect to scalar relativistic calculations is approximately 
three times larger. Therefore, it may be conjectured from these results that the large energy values 
associated with the structural perturbations may be of electronic in origin, rather than due to steric 
interaction, alone. Note that the energy associated with steric displacements of cations in oxide 
perovskites are usually of the order of few meV ( ~ 10 - 500 K). 

3.2. Electronic structure 

3.2.1. Scalar relativistic calculations To study, how the crystal field splitting (Acf) of the t 2 g - 
Cg bands evolve as a function of (f), we show in Figl2l the density of states (DOS) of SrIrOa, for 8° 
< (f> < 20°. For (j) < 4:°, the t 2 g and Cg states are strongly hybridized and results in no electronic 
gap between them. A crystal field induced gap emerges for ^ > 4°, the magnitude of which linearly 
increases with (j), as shown in the Figl^Ka). The /S.cf gap resides in the unoccupied part of the 
DOS spectra. The integrated DOS from Ep to the top of the Og band is estimated as ~ le“ per 
formula unit, consistent with the configuration of Ir ions in SrIrOa. 

Figl2 also reveals significant t 2 g band narrowing with increasing (j). The bottom of the Ir 
5d valence band at ~ —9 eV below Ep for (p = 0° linearly shifts to lower binding energies with 
increasing p, and is estimated to be at —7.5 eV for p ~ 20°. Nevertheless, the top of the t 2 g band at 
~ 0.6 eV above Ep appears to be more or less pinned, irrespective of p, for p > 4°. The variation in 
the t2g bandwidth, W(t2g), as a function of (p is shown in Figl3l(b). For the equilibrium structure, 
the t 2 g — Cg crystal field splitting was estimated as ~ 0.55 eV. Thus, IrOe octahedra tilting not only 
induce a t 2 g — eg crystal field splitting in SrIrOa, but also lead to a significant t 2 g band narrowing 
of the Ir 5d manifold. 

The atom resolved partial DOS of SrIrOa {not shown) also reveal a wide spread distribution 
of the Ir 5d and O 2p states over the entire energy region, suggesting a covalent nature of chemical 
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Figure 2. Color online: The GGA derived DOS of SrIrOs calculated with various IrOe tilt 
angles, as indicated by the legends. The inset shows the blow-up of the states near Ep and the 

" ’ ’ ’ line (black) through energy zero represents the 



Figure 3. Color online: (a) The change in the magnitude of the t2g - eg crystal field splitting of 
the Ir 5d manifold and (b) the decreasing t2g band width as a function of the IrOe tilt angle, 0 
in SrIrOa, calculated using the GGA formalism of electronic structure method. 


bonding in SrIrOa. Both bonding and antibonding states are evident in the valence band spectra, 
which are separated by a pseudo-gap like feature at —2.3 eV below Ep, i.e., the states in the energy 
range —8 < E{eY) < —2.3 represent the bonding states, while those above —2.3 eV constitute the 
antibonding states. One may note from Figl5] that the relative energy position of the bonding - 
antibonding states in SrIrOa and all non-equilibrium structures appear to be fixed, irrespective of 
the IrOe octahedral tilt angle. 

The GGA DOS find Ep on the higher side of the Og band, with ~ 4.25 St./eV/f.u., indicating 
SrIrOa to be a good metal. However, this is in sharp contrast with the experiments which 
characterize SrIrOa as a semi-metal. Furthermore, within the realms of Stoner theory of itinerant 
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Figure 4. Color online : The GGA+SOC derived DOS of SrIrOa as a function of IrOe tilt angle 
0, the magnitude of which are indicated by the legends. The inset shows the blow-up of the states 
near Ep. The vertical broken line through energy zero represents the Fermi energy. 


magnetism such high DOS at Fermi energy, N{Ep), infers to a magnetic instability. The Stoner 
product, N(Ef ) x I, where I is the Stoner factor is estimated as 1.2 in SrIrOa, with /(Ir) ~ 
0.574 [38|. Following these predictions from the GGA calculations, we carried out spin polarized 
calculations. A ferromagnetic (FM) solution was found with Ir local magnetic moment being 0.13 
Hb, which was largely associated with the splitting of the Og anti-bonding states. Beyond, the 
GGA total energy difference between the non-magnetic and ferromagnetic states was found to be 
nearly degenerate. The prediction of magnetic ordering in SrIrOa, following GGA calculations 
is in contradiction with the experiments which find no long range magnetic ordering in SrlrOa- 
Furthermore, we also extended our calculations to check for any stable antiferromagnetic solution 
with initial Ir spin configuration assigned in the antiferromagnetic A-type, C-type and G-type 
structures. However, all these calculations converged to a non-magnetic solution. These results, 
therefore, indicate that the scalar relativistic Hamiltonian although accounts for the crystal field 
gap and t 2 g band narrowing in SrIrOa, however, fails to account for the distribution of electronic 
states, in particular, that of the near Ep anti-bonding states. 

3.2.2. Effects of spin-orbit coupling Given that SOC becomes increasingly important for 5d 
elements, we show in FiglU the DOS computed by means of the GGA+SOC calculations. The 
overall features of the DOS spectra in the bonding region looks similar to that of the scalar 
relativistic calculations, while the antibonding states appear to be strongly perturbed. Comparison 
of Figl2]and Figl4]reveals that the energy distribution of the states in the vicinity of Ep, i.e., —0.5 
< E < +0.5, undergo a drastic redistribution with a valley like feature emerging with its depth 
increasing concomitantly with (f>. For the equilibrium structure {(j) ~ 16°) the N{Ep) reduces to 
1.88 St./eV/f.u, which is significantly smaller in comparison with the scalar relativistic result. The 
reduction in N(Ep) is primarily due to the partial splitting of the Jeff states of the Ir < 2 g manifold. 
According to the Jeff model of iridates, SOC splits the Ir t 2 g states into Jeff — \ and Jeff = | 
states. If one associates the Ir formal chemical valence in SrIrOs as +4, simple band filling assuming 
an ionic picture, would suggest that the low energy Jeff = | states which are a linear combination 
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of = ^ and |TOj|=| states would occupy two electrons each leaving Jeff = ^ doublet singly 
occupied. The integrated DOS from Ep to the top of the t 2 g band in SrIrOa yields ~ 1 e“per 
formula unit. However, quite different from the case of insulating Sr 2 lr 04 and Sr 3 lr 207 , the Jeff 
states in SrIrOa are found to be a superposition of both Jeff = ^ and Jef / = | states yielding 
a semi-metallic ground state. We note that the hybridization of Jeff = ^ and Jeff = f states 
in SrIrOa, more or less follows from the three dimensionality of its structure, due to which both 
in-plane and out-of-plane O 2p orbitals have considerable overlap with the Ir 5d orbitals. Note that 
in Sr 2 lr 04 and Sralr 207 , there is no orbital overlap of the Ir 5d orbitals with the out-of-plane O 
ions, primarily due to their quasi two dimensional structure. 

We also find a strong dependence of N{Ep) with (j). For example, the N{Ep) reduces to 1.06 
St./eV/f.u for (p = 20° while, for reduced tilting (p = 12° it increases to 2.39 St./eV/f.u. Apart from 
the decrease in N{Ep) with increasing <p, we also find that the energy derivative of DOS at Ep 
also increases. The latter implies to an increase in the effective mass or in other words, a decrease 
in carrier mobility (/i). Overall, the GGA+SOC calculations which find a decrease, both in N{Ep) 
and p with structural perturbations clearly suggest that the transport properties of SrIrOa can be 
controlled by tilting the octahedra, which in turn can be accomplished by strain. This is consistent 
with the strain dependent studies of orthorhombic SrIrOa films grown on different substrates. 
Tunable semi-metallic and metal-insulator transitions have been observed in experiments due to 
lattice mismatch of the orthorhombic SrIrOa with that of the substrates [30|, l31|. It is inferred 
from the experiments that in-plane compressive strain decreases the Ir 5d bandwidth. With the 
underlying symmetry of the lattice being unchanged, it can therefore be well anticipated that the 
increase in the energy of the system due to strain would be minimized mainly via the displacement 
of O ions, implying IrOe octahedral tilting. 

Moreover, in the context that octahedral tilting in SrIrOa leads to the Og-Sg band splitting 
and narrowing of the bands, then conversely, SOC via structural modifications enhance the extent 
of Ir 5d and O 2p orbital overlap. The argument follows from the fact that GGA when finds the 
equilibrium (p as 16.5°, GGA+SOG estimates it to 15.2°. Therefore, enhanced hybridization of Ir 
5d - O 2p orbitals would reflect to an enhanced dispersion of these hybridized bands. Based on 
the hybridization and Goulomb correlations analogy, our argument is consistent with the fact that 
increased SOG enhance the critical Hubbard interaction strength to drive an insulating ground 
state in SrIrOa 

3.2.3. Role of carrier concentration Having found that the structure and spin-orbit coupling play 
an important role in determining the electronic structure of SrIrOa, we investigate the role of carrier 
concentration in SrIrOa. The blow-up of near Ep region, as shown in FigHl shows that the position 
of Ep is not at the bottom of the valley, but resides along its negative slope. Therefore, based on the 
rigid band model of electronic structure it could be anticipated that electron doping, which can be 
accomplished via O defects, would slide the Ep down the valley, thereby decreasing N{Ep) further. 

In general, a low N{Ep) also indicate to higher phase stability. Therefore, such O defects modified 
electronic structure appears quite favorable in SrIrOa, which along with structural perturbations 
would increase the propensity of the system to exhibit metal-insulator transitions. 

The carrier concentration in SrIrOa is estimated to be ^ 10^° e/cmf [3|, the magnitude of 
which lies in the vicinity of the semiconductor — metal phase boundary of the solid state. In the 
rigid band model, accommodating 10^° ejcrr? (= 0.0245 eju.c.) in the DOS spectra necessitate 
a shift of Ep by 23 meV. Such an effect significantly reduce N{Ep) to 0.82 St./eV/u.c., thereby 
increasing the propensity of a metal-insulator transition in SrIrOa. Moreover, for an enhanced IrOe 
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Figure 5. Color online : The calculated DOS using the GGA+Ue// scheme (a) without and, 
(b) with SOG, for the equilibrium structure ; 0 ~ 16°. The vertical broken line through energy 
zero represents the Fermi energy. 


tilting, say (p = 20°, a carrier concentration corresponding to 10^° ejcrv? yields a of N{Ep) 0.30 
St./eV/u.c. Thus, we find that cooperative effects associated with octahedral tilting, SOC and 
carrier concentration in SrIrOa play an important role in controlling the transport properties of the 
system. 

3.2.4- Effects of Coulomb correlations We also have studied the effect of Coulomb correlations 
on the electronic structure of SrIrOa. In Fig[5Ua), we compare the GGA and GGA+Ue// DOS. 
Evidently, the most notable effect with Coulomb correlations in SrIrOa is the increase in the 
magnitude of the crystal field gap between the t 2 g and Cg bands. Also, while a small shift on 
the energy scale is also observed for the bonding states, the distribution of electronic states in 
the antibonding states remains more or less unchanged. On the other hand, the GGA+Ue//+SOC 
calculations, the results of which are shown in FigO^b), reveal no significant changes in the electronic 
structure arising due to Coulomb correlations. Thus, on a comparative length scale (0.5 - 2 
eV) between the crystal field splitting, SOC and Coulomb correlations, the effect on Coulomb 
correlations on the electronic structure appears to be very marginal. 

Moreover, to study how Coulomb correlations and SOC compete in SrIrOa, we carried out spin 
polarized calculations in the GGA+Ue// scheme for 0 < Ue// (eV) < 4 eV. In these calculations, 
an Ir local magnetic moment was computed in the FM and AFM-A and AFM-G structures, for all 
values of Ug//. For the AFM-C type ordering, an Ir local moment appeared only when Ue// >3.5 
eV. Note that the emergence of magnetic moment at the Ir sites with Ue// can be attributed to the 
increased localization of the 5d states. However, with inclusion of SOG term in the Hamiltonian, 
all calculations with Ue// < 3.5 eV, irrespective of the underlying magnetic ordering, converged 
to nonmagnetic solution. Since the origin of local moments in the GGA+Ue/ / are associated with 
localization of the 5d states, its annihilation with SOC for Ue// <3.5 eV is to be attributed to the 
delocalization (or enhanced hybridization) of the corresponding states with O 2p orbitals. However, 
we note that a stable magnetic moment (~ 0.5 fis) was found for all magnetic structures in the 
GGA+Ue//+SOC scheme of calculations when Ue// > 3.5 eV, i.e., when the magnitude of Ue// 
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exceeds the Ir 5d bandwidth. Based on these calculations, we argue that on a comparable length 
scale of interactions, Coulomb correlations and SOC compete with each other, with SOC partly 
nullifying the Ir 5d localization effects by enhancing the inter-atomic Ir and O orbital hybridization. 

4. Summary and conclusion 

In summary, we find a cooperative effect on the electronic structure of SrIrOa due to an interplay 
of its underlying lattice and SOC. Octahedral tilting of the IrOe motif in SrIrOa appears to be 
primary source of crystal field splitting of the t 2 g - Cg bands. The magnitude of the splitting scales 
with tilt angle, which also results in band narrowing. On the other hand, SOC induces a partial 
splitting of the Og bands into Je// states which reduce N(Ep), driving SrIrOa to a semi-metallic 
ground state. Taking into consideration that the tilt angle can be tuned by strain and that electron 
doping can be accomplished via controlled O defects, the propensity to drive in a metal-insulator 
transition in SrIrOa appears consistent with the experiments. Besides, on a comparative length 
scale of interactions, we find that the effect of Coulomb correlations on the electronic structure of 
SrIrOa is very nominal. Interestingly, the SOC enhanced orbital hybridization also substantiates the 
preposition that the magnitude of the critical Hubbard interaction strength to drive an insulating 
ground state in SrIrOa, increases with increasing SOC strength. 
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